Dispersive forces on bodies and atoms: a unified approach 
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A unified approach to the calculation of dispersive forces on ground-state bodies and atoms is 
given. It is based on the ground-state Lorentz force density acting on the charge and current densities 
attributed to the polarization and magnetization in linearly, locally, and causally responding media. 
The theory is applied to dielectric macro- and micro-objects, including single atoms. Existing 
formulas valid for weakly polarizable matter are generalized to allow also for strongly polarizable 
matter. In particular when micro-objects can be regarded as single atoms, well-known formulas for 
the Casimir-Polder force on atoms and the van der Waals interaction between atoms are recovered. 
It is shown that the force acting on medium atoms — in contrast to isolated atoms — is in general 
screened by the other medium atoms. 

PACS numbers: 12.20.-m, 32.80 .Lg, 42.50.-p, 42.50.Vk, 42.50.Nn 



I. INTRODUCTION 

It is well known that polarizable particles and macro- 
scopic bodies — matter whose electromagnetic properties 
are described terms of macroscopic quantities — are sub- 
ject to forces in the presence of electromagnetic fields, 
even if the (fluctuating) fields vanish on average and the 
bodies do not carry any excess charges and are unpo- 
larized. In particular, this is also the case when the 
field-matter system can be assumed to be in its ground 
state, where only quantum fluctuations are responsible 
for the forces. In this case it is common to speak of vac- 
uum forces or dispersive forces, which obviously represent 
a genuine quantum effect. Basically, it can be distin- 
guished between three kinds of dispersive forces, namely 
van der Waals (vdW), Casimir-Polder (CP), and Casimir 
forces, depending on whether forces between individual 
particles, between particles and macroscopic bodies or 
between macroscopic bodies are respectively considered 
(see, e.g., Ref. [IJ for an overview). 

Since both CP forces on atoms and Casimir forces on 
macroscopic bodies may be regarded as being macro- 
scopic manifestations of microscopic van der Waals 
forces, intimate relations between them can be expected. 
Nevertheless, quite different theoretical approaches to the 
two kinds of forces have been developed. Compared to 
the large body of work in this field, only little attention 
has been paid to their common origin and consequential 
relations between them (see, e.g., Refs. H S IE 
Moreover, the studies have been based on specific ge- 
ometries such as simple planar structures, and weakly 
polarizable matter has been considered. 

More attention has been paid to the relations between 
Casimir forces and van der Waals forces, but again for 
specific geometries and weakly polarizable matter (see, 
e.g., Refs. @, 0, II S ES El 03). Quite recently, a 
very general relation between CP forces (as calculated 
within the frame of macroscopic QED) and multi-atom 
van der Waals forces has been established [l3|, where it 
has been shown that the CP force acting on an atom in 
the presence of a dielectric body of given permittivity can 



be regarded as being the sum of all many-atom van der 
Waals forces with respect to the atoms of the body. 

In this paper we develop, within the framework of 
QED in linearly, locally, and causally responding magne- 
todielectric media, a unified approach to the calculation 
of dispersive forces acting on ground-state macro- and 
micro-objects. Since the origin of any electromagnetic 
force is the Lorentz force acting on appropriate charges 
and currents, we first consider the ground-state expec- 
tation value of the Lorentz force density acting on the 
charge and current densities attributed to the polariza- 
tion and magnetization fields of linear media, taking fully 
into account the noise polarization and noise magneti- 
zation that are associated with absorption. From the 
ground-state Lorentz force density obtained in this way, 
the force acting on an arbitrary body or an arbitrary part 
of it can be then obtained by integration over the respec- 
tive volume. Applying the theory to dielectric systems, 
we present a very general force formula, the applicability 
of which ranges from dielectric macro-objects to micro- 
objects, also including single atoms, without restriction 
to weakly dielectric material. In particular, this formula 
enables us to extend the well-known CP-type formula for 
the force acting on a weakly dielectric (micro-) object to 
an arbitrary one. 

The paper is organized as follows. In Sec. [H] the the- 
ory of CP forces acting on single ground-state atoms in 
the presence of magnetodielectric bodies is recapitulated. 
The ground-state expectation value of the Lorentz force 
density in linear magnetodielectric media is calculated in 
Sec. IIIII On this basis, in Sec. HVI the force acting on an 
arbitrary dielectric body or a part of it is derived. The 
theory is applied in Sec.[V]to the calculation of the force 
acting on a dielectric micro-object, the limiting case of 
a single atom (which can be either an isolated atom or 
a medium atom) is considered, and contact to earlier re- 
sults found for planar structures is made. In Sec. IVII it 
is shown that the theory can also be used to study the 
van der Waals interaction between two atoms. Finally, 
a summary and some concluding remarks are given in 
Sec. I VIII For the sake of clarity, some of the derivations 
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are given in appendices. 



II. CASIMIR-POLDER FORCE 

The CP force acting on a ground-state atom in the 
vicinity of arbitrary magnetodielectric bodies that lin- 
early, locally, and causally respond to the electromag- 
netic field can be regarded as being a conservative force. 
Hence it can be given by the negative gradient of a po- 
tential which in the lea ding order of perturbation theory 
reads (see, e.g., 0HIIM1G1) 

P (at) (r) = _^0 r d ^2 a( . c)VTrG (S) (r)r . 0j (1) 
27T Jo 

where r is the position of the atom, a(z£) is its polariz- 
ability, and G^ s ^ (r, r', i£) is the scattering part of the clas- 
sical retarded Green tensor G(r, r', i£) taken at imaginary 
frequencies. Note that the scattering part of the Green 
tensor contains all the necessary information about the 
configuration of the magnetodielectric bodies, whose elec- 
tromagnetic properties are characterized by the electric 
and magnetic susceptibilities that are complex functions 
of frequency and may vary in space. The (translationally 
invariant) bulk part of the Green tensor, which would di- 
verge in the coincidence limit, is not needed in Eq. l[T]l. 
because it cannot contribute to the force. Equation 0J, 
which can be derived 0] on the basis of exact quantiza- 
tion of the macroscopic electromagnetic field in linearly, 
locally, and causally responding media strictly ap- 
plies to isolated atoms, but not to medium atoms nor 
to guest atoms in a substrate medium. Note that the 
atomic ground-state polarizability a(u>) in leading-order 
perturbation theory, 

k k 

features poles on the real frequency axis due to the ne- 
glect of level broadening. If necessary, the correct re- 
sponse function properties pol | may be restored by means 
of an appropriate limit prescription, viz. 

a{uj) ~ lim Y — "* . . (3) 

In order to apply Eq. Q to a dielectric body of volume 
Vm, let us consider instead of a single atom a collection of 
atoms that are (strictly) contained inside a space region 
of volume Vm, and let us add up the individual forces 
as given by Eq. Q. Since the mutual interaction of the 
atoms is completely ignored in this way, it is clear that 
this method gives only the lowest-order approximation 
to the total force. If the number density of the atoms 
(defined on a suitably chosen macroscopic length scale) 
is denoted by rj(r), the total force in this approximation 



reads 

F = Jjfo f d 3 r r d ^ (r)a(lC)V TrG< s >(r,r 
m Jv M Jo 

(4) 

Since the validity of Eq. Q obviously requires sufficiently 
weakly polarizable atoms and/or a sufficiently low num- 
ber density of atoms, the collection of atoms can be 
viewed as dielectric matter of volume Vm and small sus- 
ceptibility 

»(r, i£) = r)(r)a(i£)/E , (5) 

which implies that the permittivity of the overall system 
has slightly been changed by 6s(r, i£) = XM.( r ,iQ- In 
particular, applying Eq. @| to a dielectric micro-object 
whose number density of atoms is constant over the small 
volume Vm, we obtain the force 

F = - Vuv ^ J°° d£ ea (i0 VTr (r, r, if). (6) 

Clearly, application of Eq. Q to dielectric bodies (in- 
cluding micro-objects) that are dense and/or consist of 
strongly polarizable atoms becomes questionable. Hence, 
the problem of a generalization of Eq. Q to arbitrary 
dielectric bodies or parts of them arises. As we will see, 
a satisfactory answer can be given on the basis of the 
ground-state Lorentz force density in media. 

III. GROUND-STATE LORENTZ FORCE 

As shown in Ref. [2l|, the Casimir force between (lin- 
early, locally and causally responding) macroscopic bod- 
ies can be regarded as the expectation value of the 
Lorentz force acting on the charges and currents at- 
tributed to the polarization and magnetization of the 
bodies. To be more specific, the charge and current den- 
sities p(r) and j(r), respectively, which the Lorentz force 
density acts on are given by 

/>oo 

p(r) = / dojp(r,Lu) + H. c. (7) 
Jo 

and j (r) accordingly, with 

£(r,w) = -£ V-{[e(r,w)-I]E(r,w)} + (^)- 1 V-j N (r,w) 

(8) 

and [k = M ] 

j(r,w) = -iue Q [e(r,u) - l]E(r,w) 

+ V x {/s [l-K(r J w)]B(r,w)}+j N (r,w). (9) 

Here, j (r, u;) is the (fluctuating) noise current density 
that acts as a Langevin noise source in the (macro- 
scopic) Maxwell equations in the domain, e(r,u) and 
K _1 (r, u)) are respectively the (complex) permittivity and 
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permeability, and E(r, u>) and B(r, u>) are the (positive) 
frequency components of the medium-assisted electro- 
magnetic field: 



E(r, w) = ifiooJ I dV G(r, r', u) ■ j N (r', u), (10) 
B(r, W )= At0 Vx /dVG(r,r>)-L(r>) (11) 



(for details of the quantization scheme, see Ref. Il9j). 
The expectation value of the Lorentz force, 



d 3 



/■ / dw / du/(£(r,w)Et(r',u/) 

JO 



+ j(r, W )xB'(rV'))r 



(12) 



taken with respect to the ground state of the linearly 
interacting field-matter system then yields the Casimir 
force acting in the zero-temperature limit on the material 
in a chosen spatial region of volume Vm ■ Again, divergent 
bulk contributions must be discarded in the limit r'— >r, 
as they would correspond to an unphysical "self-forces" 
of the respective volume elements. 

The Casimir force as given by Eq. (|12fl [together with 
Eqs. (|5|)- (|ll|l ] can be equivalently rewritten as a surface 
integral over a stress tensor, 



da • T(r) 



<3V M 



where (at zero temperature) 

T(r) = lim [S(r,r') - ±lTrS(r,r') 

r'— >r 

together with 



(13) 



(14) 



S(r.r') 



d£ 



G (S) (r,r',z6 



+V x G (s) (r,r',^) x V' 



(15) 



Here and in the following, 1 denotes the unit tensor and 

xV' is meant to act to its left. With regard to the con- 
troversial view held in Ref. [22] , it should be emphasized 
that Eq. (JT3J together with Eqs. and JTJjJ yields the 
genuine electromagnetic Casimir force acting on bodies 
or pieces of them. Clearly, in mechanical equilibrium the 
Casimir force is balanced by additional internal or ex- 
ternal (mechanical) forces that are not included in the 
equations considered here. 

Let us return to Eq. (TJJ together with Eqs. (|5|)- ljTT)l . 
It is not difficult to show (Appendix^) that Eqs. © and 
© can be rewritten as 



p{r,u) = -V- / dVG(r,r».L(r' )W ), 



(16) 



j(r,w) = V x V x 



dVG(r,r».j N (r', W ), 
(17) 



and that the relation 
(j N (r,w)j^(r',w')) 



8{w-uj')\ %rlm e(r,Lu)16(r -r') 

Mot L tr 



- V x [Im «(r, w)15(r - r')] x V \ (18) 



holds for the ground-state expectation value. Employ- 
ing standard properties of the Green tensor (see, e.g., 
Ref. 0), from Eqs. (EH), OH and jI fy -fT5 )l it follows 
that 

(p(r,u)&(r',u>')) = -^( W - W ')VdmG(r,r» 

— 7T Cr 

(19) 

and 

(|(r, W )Bt(r',u/)) 
= —S(u} - V x V x ~]ImG(r,r» x V, 

7T \ C 2 J 

(20) 

which can be used to express the integrand of the volume 
integral in Eq. (|12|l — that is, the Casimir force density — 
in terms of the (scattering part of the) Green tensor 
solely. Extension of the theory to include thermal states 
is straightforward. 



IV. FORCE ON DIELECTRIC BODIES 

To illustrate the theory, we first apply Eq. (|12|l [to- 
gether with Eqs. (|19|) and (|20|) ] to the calculation of 
the Casimir force acting on a dielectric body in some 
space region of volume V . We want to study — mostly in 
parallel — the two cases sketched in Figs. Handle namely 
(i) an isolated body (Fig.[I}, an d (ii) a body that is an in- 
ner part of some larger body (Fig. [2J ■ In both cases, arbi- 
trary magnetodielectric bodies are allowed to be present 
in the outer region Vb in the figures. 

Let e(r, w) be the permittivity of the system in the ab- 
sence of the dielectric matter in V and assume that it 
changes to e(r, w) + Ae(r, ui) when the additional dielec- 
tric matter is introduced into the initially empty space 
region V. Here and in the following, AA denotes the ex- 
act change of a quantity A[e(r, w)] due to a given (not 
necessarily small) change of the permittivity, Ae(r,w), 
whereas the notation 5 A is used to indicate the familiar 
first-order variation of A produced by a small variation 
6e(r,Lj). Further, let Gy(r,r',u;) and G(r,r',u;) be the 
Green tensors of the system in the cases where the di- 
electric matter inside the space region V is present and 
absent, respectively, with both of them taking into ac- 
count the magnetodielectric bodies in the space region 
Vb in Figs. ^ and [21 We may then write 



Gy(r,r» = G(r,r» + AG(r,r', 



(21) 
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FIG. 1: A dielectric body M of volume Vm inside an empty- 
space region of total volume V. There may be arbitrary mag- 
netodielectric bodies in the outer region Vb. 




FIG. 2: A dielectric body M of volume Vm that is an inner 
part of a larger dielectric body of volume V. There may be 
arbitrary magnetodielectric bodies in the outer region Vb. 



where the change of the Green tensor, AG(r,r',w), due 
to the introduction of the dielectric matter into V obeys 
the Dyson-type equation 

uj 2 f 

AG(r,r',w) = - / d 3 s G(r, s, uj) ■ Ae(s, w)Gy(s, r', w), 

(22) 

with the integral running over the space region V (at 
most). 

To calculate the changes of the expectation val- 
ues (/5(r,w)Et(r',w')} [Eq. £2)] and (j(r, w)Bt(r',w / )) 
[Eq. pb)l] (r G V), we note that for r G V the Green ten- 
sor G(r,r',w) satisfies the same differential equation as 
the free-space Green tensor, i.e., 

(v x V x ~} G(r, r', u) = 16(r - r') (r G V), 
from which it follows that 

UJ 2 

- T V-G(r,r',u) = -V5(r-r') (r£V). (24) 



Hence from Eqs. (|22() and l|24|l we derive 

V • ImAG(r,r',w) 

= -V-Im [Ae(r,w)G v (ry,w)] (r G V), (25) 

and thus from Eq. (|T§)) the change A(p(r, w)E t (r', w')) 
is found to be 

A(p(r, W )Et(r',c')) = --^^-c') 

X V-Im [Ae(r,w)Gy(r,r',cj)] (r € V). (26) 

The change A(j(r, cj)Bt(r', a/)) (r G V) can be found in 
a similar way. Making use of Eqs. 1221) and l|23|) . from 
Eq. lj2U|) we obtain 

A(j(r,o;)Bt(r',a; / )) = --^^- W ') 

xIm[AE(r,w)Gv(r,r»]xV' (r G V), (27) 
which implies that 

A(j(r,c) xBt(r>')) 

= - ^*(w - c/)Im{A£(r,w)V'Tr[Gy(r,r',a;)l 
-Ae(r,w)V'-Gy(r,r',w)} (r G V). (28) 

Using Eqs. (|26|l and (J2SJ, we can now easily calculate, 
according to Eq. (jl2JI . the Casimir force acting on a di- 
electric body of volume Vm and dielectric susceptibility 
XM(r, oS) = Ae(r,w) (cf. Figs. [l] and 0) as 




x Im / d'V {xM(r,w)V'Tr[Gy(r,r',w)] 
- (V + V) • XM (r,w)G v (r,r»} r ,^. (29) 

The reciprocity property of Gy(r, r',u/) implies that 
Tr Gy (r, r', w) is symmetric with respect to r and r'. 
Thus we may rewrite Eq. I|29|) as 

F = ^ / dww 2 

x (im / d 3 r X M(r,w)VTr[Gy(r,r', W )] r ^ r 

-21m / da-xu(r,u})[G v (r,r',cj)] r '^ r \. 
Jav M ) 

(30) 

Further, on recalling the analytic properties of the in- 
tegrands as functions of (complex) uj, we may employ 
contour integral techniques to represent Eq. Q3U|1 in the 
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form of 

X (/ d 3 r XM (r,^)VTr[Gy(r,r',iO]r^r 

-2/ da-XM(r,»0[Gv(ry,»OW}. (31) 

In Eqs. (|29() - (|31|l the coincidence limit r'— >r has again to 
be performed in such a way that unphysical "self-force" 
contributions are removed. 

If the body consists of homogeneous dielectric matter, 
one can simply replace the Green tensor [Gy (r, r', uj)] r /^ r 

with its scattering part G v (r, r, w). In the case of in- 
homogeneous matter this replacement should be done 
point-wise. That is to say at each space point r, the 
Green tensor for the corresponding bulk material must 
be subtracted from Gy(r, r', to) in the limit r' — >r. In the 
case of an isolated body (cf. Fig. 2} the surface integral 
in Eqs. I|3U|) and l|31fl can be dropped, hence 




x f d 3 rxM(r,iOVTr[Gv(r,r',z£)] r ^ r . (32) 
Jv M 

Clearly, the surface integral must not be dropped in the 
case where the body is an inner part of a larger dielectric 
body (cf. Fig. . 

Equation i|32|) is the desired generalization of Eq. J3J. 
In contrast to Eq. 0}, it represents the exact force acting 
on a dielectric body of given permittivity, since xm(f, «£) 
is not restricted to small values anymore. Correspond- 
ingly, the Green tensor in Eq. (1321) is the one that takes 
the presence of the dielectric body fully into account, 
whereas in the Green tensor in Eq. 0} the presence of 
the dielectric body is not considered. Hence in contrast 
to Eq. Q}, Eq. (|32|l includes in the calculation of the 
Casimir force that acts on a dielectric body the body's 
retroaction on the electromagnetic ground-state noise of 
the residual system. 

To make this explicit, one can expand the full Green 
tensor Gy(r, r',i£) in powers of Ae(r, i£) by using the 
iterative solution to Eq. (|2~2"|) [together with Eq. JH}]. 
Inserting the resulting Born series for Gy(r, r',i£) in 
Eq. (|31|l . one obtains the corresponding expansion of the 
Casimir force F in powers of Ae(r, In particular, 
truncating this expansion at the term linear in Ae(r, i£) 
[Ae(r, i£) 1— > Se(r, ig)], i.e., replacing Gy(r, r', i£) with its 
zeroth-order approximation G(r, r', we simply obtain 

x ( / d 3 r XM (r, iQVTr [G(r, r', 

-2/ da. XM (r,^)[G(r,r',iO]r^r). (33) 

JdVm > 



In this case, the prescription for taking the coinci- 
dence limit of the Green tensor simply consists in 
the replacement [G(r, r', i£)] r '-»r | — * G^ s ^(r, r, where 
G( s )(r, r', i£) is the scattering part of the Green tensor 
G(r, r',i£) in the absence of the dielectric matter in V 
(cf. Figs. n]aud[5}. Note that if the surface integral can 
be dropped (i.e., if the case sketched in Fig. ^ is con- 
sidered), then Eq. I|33[) becomes identical with Eq. 0}. 
It should be mentioned that inclusion in Eq. (|33() of the 
higher-order terms of the Born series of Gy (r, r',i£) gen- 
erates an increasing number of many-body corrections. 

It may be also informative to examine the effect of the 
change Ae(r, uj) not only on the level of the expectation 
values in Eqs. Ijl9(l and J20}, but more directly on the 
level of the corresponding operators. This is outlined in 
Appendix for weakly dielectric material [Ae(r,u>) 
6e(r,u))]. In this context, an alternative derivation of 
Eq. l|3l^l is given. 

In view of Eq. (|13|l . Eq. (|31|l can be also given in the 
form of a surface integral, where the stress tensor T(r) 
can be found by studying [in a similar way as in the 
derivation of Eq. 1|31|) ] the change induced in Eqs. Ijl4(l 
and (|15|) by the change Ae(r, More easily, we may di- 
rectly derive it (up to irrelevant transverse contributions) 
from Eq. ll3"T|) as 

T(r) = — / d££ 2 XM(r,i£)[Gy(r )M £)] r ^ r 
nc J 

2ttc 2 J a J Vm 

x XM(r',iO , y 
47r r — r 

which makes obvious the fact that the stress tensor de- 
pends on the permittivity in a spatially non-local way in 
general. With respect to the volume integral in Eq. I|34|l . 
r should be thought of as being infinitesimally outside 
Vm if necessary. The formulation of the force in terms of 
the stress tensor is particularly advantageous in the case 
of homogeneous material. From Eq. (|31l) it is easily seen 
that for a homogeneous body that is an inner part of a 
larger body (cf. Fig. [2} the somewhat cumbersome stress 
tensor (|34f> can be replaced with 

h f°° 

T(r) = ^ / d££ 2 XM^){ilTr[Gy(r,r'^)] 

nc Jo 

-Gy(r,r'^)} r ,_ r , (35) 

and for an isolated homogeneous body (cf. Fig. ^ it 
follows from Eq. Q32JI that it can be replaced with 

T(r) = -tt^ dee 2 XM(^)lTr[Gy(r,r',iO] r ^ r . 

(36) 

Furthermore, the assumed homogeneity then implies that 
we may let [Gy (r, r',i£)] r /_> r 1— > Gy ^(r, r, i£) in Eqs. l|33|) 
and (|36|l . Needless to say that replacing the full Green 
tensor Gy(r, r',i£) as appearing in Eqs. I|34 |l -ip5 )l with 



6 



the zeroth-order approximation G(r,r',i£) yields again 
the Casimir force in the case of weakly dielectric material. 



V. FORCE ON MICRO-OBJECTS AND ATOMS 

Since nothing has been said about the spatial exten- 
sion of the bodies under consideration, the applicability 
of Eqs. (|31|) and IMl'l) ranges from dielectric macro-objects 
to micro-objects, even including single atoms. Let us 
consider a dielectric body that may be thought of as 
consisting of distinguishable (electrically neutral but po- 
larizable) micro-constituents frequently called atoms or 
molecules within the framework of molecular optics. We 
may then assume the validity of the Clausius-Mosotti 
relation Ell^l. 



XmO, u) = e 1 r?(r)a(w)[l - r)(r)a(uj) / (3e )\ 1 

= eg Vr)aM [1 + Xm (r, w)/3] , (37) 

where a(u>) is the polarizability of a single micro- 
constituent and ?y(r) the number density of the micro- 
constituents (referred to as atoms in the following). It is 
worth noting that there is no need here — in contrast to 
Eq. 0} — to regard a{ui) as being calculated in the lowest 
(non-vanishing) order of perturbation theory according 
to Eq. J5J). It can be shown (Appendix that Eq. 1 (57)1 
is consistent with the requirement that both a(u>) and 
XmO*, u>) be Fourier transforms of response functions iff 



77(r)a(0)/(3£ ) < 1. 



A. Isolated micro-object 



(38) 



Let Vm be the small volume of an isolated dielectric 
micro-object (cf. Fig. which a dielectric susceptibil- 
ity xm(w) of Clausius-Mosotti- type can be ascribed to. 
Combining Eq. '32fl with Eq. '37fl an d assuming that, 
due to the smallness of Vm, the scattering part of the 
Green tensor can be taken out of the space integral at the 
(appropriately chosen) position r of the micro-object, we 
derive the force acting on the micro-object to be 

F = -V MV ^ [°°dU 2 a(iO [1 + |xm(»0] 
Z7r Jo 

x VTrG^ S) (r,r,^). (39) 

Recall that in the case under study the replacement 
[Gy(r, r , j£)] r /_n. i— ► G v '(r,r,i£) can be made. Equation 
(|39|l . which generalizes Eq. JJJ, differs in two respects 
from Eq. ©. Firstly, its validity is no longer restricted 
to weakly dielectric matter. Secondly, it takes into ac- 
count the dependence of the force on the shape of the 
micro-object. 

In contrast to Eq. © , the force as given by Eq. I|39|) in- 
cludes all-order multi-atom van der Waals interactions of 
the micro-object, as may be seen by expanding the Green 



tensor Gy(r, r',i£) in powers of Xm(*£) (cf- Ref. [13j ). If 
they are disregarded, Eq. 1(55)) reduces to (Appendix [DJ 



-V M il 



— /°°d££ 2 a(^)VTrG (s) (r,r,z£), ( 40 ) 
t Jo 



which, as expected, is nothing but Eq. @ — only the term 
linear in a(i£) contributes to the force. The force in this 
limit is simply the sum of the forces acting on the atoms 
due to the presence of the external bodies (region Vb in 
Fig. P). Hence, F( at ) = (V M ?7) _1 F is the force acting on 
a single ground-state atom, that is to say, we are left 
exactly with the formula for the CP force as given by 
Eq. (JTJ, with the exception that now the atomic polar- 
izability is the exact one rather than the perturbative 
expression given in Eq. J3J). 



B. Micro-object that is an inner part of a larger 
body 



Let now Vm be the small volume of a dielectric micro- 
object that belongs to a larger body of volume V of the 
same atoms (cf. Fig- EJ) - Under assumptions analogous to 
those leading from Eq. to Eq. from Eq. (PU [to- 
gether with Eq. (J37J)] we obtain the following formula for 
the (shape-dependent) force acting on the micro-object: 

F = -Vuv ^ J™ d£ e a(iO [1 + IxuiiO] 



x V- 



ilTrG^r.r.iO-G^wO . (41) 



Equation l|41l) differs from Eq. (|39|l in the second term 
in the square brackets in the second line. This difference 
can be regarded as reflecting the fact that — in contrast to 
Eq. I|39|) — the force acting on the micro-object is screened 
by the residual part of the body. 

When the multi-atom van der Waals interactions of the 
body (of volume V) are disregarded, then, in a way quite 
similar to that outlined in the derivation of Eq. 14UI) in 
App. [D] Eq. H41|) can be shown to reduce to the term 
linear in the atomic polarizability, 



F = - VuT] 

7T 



da 2 am 



x V 



ilTrG (s) (r,r,z£) - G (s) (r,r,^)l • (42) 



Recall that in this approximation Gy^(r, r, i£) can be 
replaced with (r, r, From Eq. (|42|l it then fol- 
lows that F' at ) = (Vm^) _1 F can be regarded as the 
screened CP force acting on an atom of a weakly dielec- 
tric medium. 

To make contact with earlier results, let us apply 
Eq. H42f> to the atoms of a weakly dielectric medium (cor- 
responding to the region V in Fig.[2J) in front of a laterally 
infinitely extended magnetodielectric planar wall (corre- 
sponding to the region Vb in Fig.[2Jl, which is assumed to 
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extend from some negative z value up to z = 0. Using the 
explicit form of the Green tensor for planar multi-layer 
structures (see, e.g., [H|2i|), we may write its scatter- 
ing part for coincident spatial arguments in the (empty) 
space region z > as 

-P 2 e q e q ] +r s _k 2 e s e s }, (43) 

with k = fc(w) = w/c, q = |q|, (3 = (3(oj, q) = (k 2 - q 2 ) 1 ' 2 
and orthogonal unit vectors e q = q/q, e z = Vz, and e s 
— e q x e z . The effect of the (multi-layered) wall is de- 
scribed in terms of the generalized reflection coefficients 
r°_ =rZ{co,q) (a = s,p), which in the simplest case of an 
internally homogeneous, semi-infinite wall reduce to the 
usual Fresnel amplitudes. From Eq. I|43fl it then follows 
that 



TrG^( 




(44) 

Substitution of Eq B4l in to Eq. (|4T1 leads to the well- 
known expression 011 E [H HI for the CP force 
acting on a single ground-state atom in front of a planar 
wall. The screened force acting on a weakly dielectric 
medium atom is obtained by substituting Eqs. i|43[) and 
gl|| into Eq. gSJ. The result is (/3 = in) 

F^(z) = (V M r 1 r 1 F(z) 

= e^J o °°dU 2 aml°°d qq e- 2 «*(rL-r p _). (45) 

It fully agrees with the result found by calculating the 
Casimir stress (|14fl [together with Eq. (|15f) ] in a dielectric 
layer of a planar multi-layer structure and performing 
therein the limit to weakly dielectric matter jj, 0] . 

VI. VAN DER WAALS INTERACTION 
BETWEEN TWO ATOMS 

Equation (|3*T|) can also be regarded as a basic equa- 
tion for calculating the force between two (ground-state) 
atoms. For this purpose, let us consider the small change 
SF of F in Eq. (|31ll due to a small change <5xi(r,i£) of 
the susceptibility xm(f, i£) and a small change #X2(r, i£) 
of the susceptibility xb (r , i£) (of one of the bodies) in the 
region Vb in Fig. J2J). In particular let us assume that 
Xm(i",«£) only changes inside the region Vm- Recalling 
Eq. (I22|) . it is not difficult to calculate 6F up to second 
order in <!>Xfc(i", ££) (k = 1, 2) and pick out the term 5i2F 
that is bilinear in <5xi(r, i£) and <!>X2( r , «£): 

h r°° r 
Si2F=—^ / d^ 4 / d 3 r6 X i(r,*0 
Z7TC Jo Jv M 

x f d 3 S( 5x2(s,iOVTr[Gy (r,s,tf).Gv(8,r,iO], 

(46) 



where the Green tensor Gv(r%, T2, i£) refers to the system 
before the susceptibilities have been changed. Note that, 
since we are dealing with the interaction between two 
well-separated space regions, the problem of removing 
"self" -force contributions does not arise here. 

Now let us suppose that the small changes 5\i (r, i£) 
and <5x2( r >*£) result from the introduction into the sys- 
tem of additional atoms, say impurity atoms, of type 1 
and type 2, respectively. The (body- assisted) force act- 
ing on a type-1 atom at position ri due to its interaction 
with a type-2 atom at position Y2 is then evidently ob- 
tained, in first order of their polarizabilities ai(i£) and 
a2(i£), from the "crossing term" <5i2F as 

fc 2 poo 

x VilirlGvCri.ra.iO-Gvfra.ri,^)], (47) 

which is in full agreement with previous calculations of 
the van der Waals interaction between two atoms IH, 
l20l | . Recall that Gy(ri, r%, i£) is the Green tensor for the 
material system that was present before the introduction 
of the additional atoms. 

Disregarding local-field corrections, one may insert in 
Eq. (|47l) the Green tensor for the unperturbed host me- 
dia. In particular, the force between two atoms embed- 
ded in a homogeneous (dielectric) background medium 
is then obtained by identifying Gy(ri,r2,i£) with the 
well-known bulk-medium Green tensor. Note that in this 
case the same formula for the force can be obtained by 
ba sing the calculations on the Minkowski stress tensor 
H<3] Choosing in Eq. |g7J| the free-space Green 
tensor, we recover the van der Waals interaction between 
two atoms in otherwise empty space. It should be pointed 
out that F^ t} and F ( 2 f } obey the lex tertia F^ t} = -F 2 f } 
if the Green tensor is translationally invariant [Gy(ri + 
v,r2 + v, i£) — Gy(ri,r2, as it is the case for the 
two atoms being in bulk material or in free space. Since 
Eq. (|47ll describes the atom-atom force in the presence of 
arbitrary macroscopic bodies, it is clear that the atomic 
positions ri and r2 are not physically equivalent in gen- 
eral. 



VII. SUMMARY AND CONCLUSIONS 

Within the framework of macroscopic QED in linearly, 
locally, and causally responding media, we have shown 
that dispersive forces acting on (ground-state) macro- 
and micro-objects — including single atoms — can be cal- 
culated in a unified way on the basis of the Lorentz force 
density that acts on the charge and current densities at- 
tributed to the polarization and magnetization of the 
media. Although the examples considered in Sees. IIVI - 
IVII refer to dielectric objects, the basic formulas given 
in Sec. IIIII can also be used to include in the calcula- 
tions magnetic properties of the matter. Inclusion in the 
theory of non-locally responding media would require an 
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extension of the underlying quantization scheme, which 
may be a subject of further studies. 

We have derived very general formulas for the force act- 
ing on a dielectric body or a part of it — formulas which 
apply to arbitrary geometries and whose validity is not 
restricted to weakly dielectric matter. For locally re- 
sponding dielectric matter that may be regarded as con- 
sisting of atoms in the broadest sense of the word, the 
permittivity can be assumed to be of Clausius-Mosotti- 
type. In this way, all relevant many-atom van der Waals 
interactions of the involved matter can be included in the 
force to be calculated. 

As already mentioned, the applicability of the theory 
ranges from macro-objects to micro-objects. Commonly, 
the force acting on a (weakly) dielectric micro-object is 
calculated in the spirit of a simple superposition of CP 
forces acting on independent atoms. The present theory 
enables one to systematically include in the calculation 
both the dependence of the force on the shape of the 
micro-object and, at the same time, the contributions 
to the force due to many-atom interactions of atoms of 
the micro-object, without restriction to weakly dielectric 
matter. 

If the micro-object reduces to a single atom, the well- 
known formula for the CP force on a single atom is re- 
covered. It is worth noting that not only the force acting 
on an isolated atom can be obtained, but also the force 
on a medium atom. For a medium atom, the CP force 
is screened due to the presence of neighboring medium 
atoms, while there is of course no such screening in the 
case of an isolated atom. Moreover, the basic formulas 
can also be used to study the body-assisted van der Waals 
interaction between atoms. 
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Eq. iTATTi . 



-^V-[e(r,ui)G(r,r',ui)] = -V6(r-r'), 



(A2) 



we easily see by straightforward calculation that 
Eqs. i|Pp|) and (fTTfl hold. According to Ref. the noise 
current expressed in terms of the noise polarization and 
the noise magnetization, 



(A3) 



can be related to bosonic fields f\(r, ui) and fUr, cj) 

(X = e,m), 

[/ Afe (r, W ),/^(rV)] = 5 kl 6 xx ,S(r-r')S(ui-ui'), (A4) 
by means of the relations 

P N (r,w) = i [fi£ Ime(r, ui)/n] 1 f e (r,w), (A5) 

1/2 " 

M N (r, ui) = [-7wcoIm«(i-, uj)/it] f m (r,co), (A6) 

Note that the f\(r, ui) and f|(r, ui) play the role of the 
basic variables of the combined system composed of the 
electromagnetic field and the (linear) medium. Insert- 
ing Eqs. I|A5(1 and l|A6(l in Eq. (|A3|> and making use of 
Eq. (|A4() we find that J N (r, ui) and jt (r, ui) obey the com- 
mutation relation 

[2 Nk( r, W )Jt i(r ', W ')]=^^--') 

,2 



— 2-\/lme(r, ui) 15{v— r')^/lme(r', u/j 



Vx 0m/c(r,w) lS(r — r')y/lmK(r',ui') x V' 



kl 

(A7) 



which immediately implies the expression for the ground- 
state expectation value (j (r, w)j (r', u/)) as given in 
Eq. dl). 



APPENDIX B: ALTERNATIVE DERIVATION OF 
EQ. (1351) 



APPENDIX A: DERIVATION OF EQS. (tT51)-(tT51) 

To express p(r,u>) and j(r, ui) as defined by Eqs. (jBJ 
and @, respectively, in terms of j (r, ui), we first insert 
Eqs. ||TJU and Jill) in Eqs. © and ©. Taking into ac- 
count that the Green tensor obeys the differential equa- 
tion 

ui 1 

V x k(t, w)V x G(r, r', ui) 5-efr, w)G(r, r', ui) — 15(r — r') 

c z 

(together with the boundary condition at infinity) as well 
as the relation that follows by taking the divergence of 



Let us consider a small variation fe(r, ui) and the cor- 
responding (first-order) changes of the operators p(r, u>), 

j(r, ui), E(r, ui), and B(r, ui), and assume that 8e{r,ui) 
is different from zero only in the volume V [cf. Fig. |2J. 
From Eqs. ©-{HI it follows that 

5pir,ui) = -eoV • {5e(r, w)E(r, ui) 
+ [e(r,w) - l]<5E(r,^)} + (iui)- x V ■ 6} N (r,ui), (Bl) 

<5j(r, ui)=- iuie {Se(r, w)E(r, ui) 

+ [e(r,ui) - l}6±(r,ui)} + 5l N (r,ui), (B2) 
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SE(r, u) = ifioui J d 3 s [G(r, s, us) ■ c5j N (s, w) 

+ 5G(r,s,w).j N ( B) w)], (B3) 

<5B(r, u) = MoV x J d 3 s [G(r, s, w) • £j N (s, w) 

+ «G(r,s J w).j N ( B) w)]. (B4) 

Since the state space attributed to the dynamical vari- 
ables fx(r, ui) and f|(r, to), in terms of which the electro- 
magnetic quantities are thought of as being expressed, 
can be regarded as being independent of the chosen per- 
mittivity (and/or permeability) |Tflj . we may apply the 
rule (8 ■ ■ ■} = 8(- ■ •} when calculating expectation values. 
Thus, combining Eqs. ©-(EJ with Eqs. l|BT) - l|B4|) . we 
derive 

5(p(r J w)E t (r',w / )> = -e V-[*e(r,a;)^(r ) «)E t (r / 1 w , )>] 



ij V ■ [(*j N (r,a;)E t (r' >W ')> + <j N (r, ^E f (r>')>] 



(B5) 



and 



«{j(r,w)B + (r / ,w')> = -tWE fc(r, w)^(r, w)B + (r / , w')> 



[<C(r, W )B t (r', W ')> + (Ur^mVV)}] 



(B6) 



where terms that vanish when r, r' G 1/ have not been 
quoted [e(r, a;) = 1 in V]. Now from Eq. i|12[l together 
with Eqs. and JEBJl and Eqs. CHI), JTTJl, d, (|B3)l . 
and (|B4|> we may easily calculate 5F. On recalling stan- 
dard properties of the Green tensor, we derive 



sf = + <y( 2 )p, 



(B7) 



where 



ft '°° <- 2 



(5 (1) F = — / duj 



2tt Jo 



I d 3 rfe(r,w)VTr[ImG(r,', 

JVm 

-2 1 da • 8e(r, w)Im [G(r, r', w)] r ^ r J. (B8) 



arises from the first terms on the right-hand sides of 
Eqs. CEH) and ljBB]l . and 



2- 



dcj ■ 



/ d 3 r [Imfe(r,cj)]VTr[G*(r,r',w)] r /. 
-2/ da- [Imfe(r,w)][G*(r,r',w)] r 



arises from the second terms on the right-hand sides 
of Eqs. (|B5|I and l|B6|l . Note that in the derivation of 
Eq. I|B9|) . the relation 

<4(r,^(r>')} = <jyr, W )5jy V)) 

= <S(w - u')% Im fefr, w)l«y(r - r') (BIO) 

2/io7T 

has been used, which follows from Eq. Q18[l. Combining 
Eqs. (E3, (E83), and (JbU yields 



£F = — /"du,^ 



2tt ./o c 2 



x Imj / d 3 rfc(r,cj)VTr[G(r,r',w)] r ^ r 

-2 [ da-fc(r,w)[G(r,r',cj)] r /^ r l. (Bll) 

Changing in Eq. 1B11|I the real-frequency integral to an 
imaginary-frequency integral in the usual way, we just 
arrive at Eq. |J53J| [6e(r,ig) i-> x M (r,?£); <SF F]. Note 
that only the real parts of Eqs. (|B8|1 and l|B9|) contribute 
to Eq. (|B11() . while the imaginary parts drop out. 



APPENDIX C: CLAUSIUS-MOSOTTI 
SUSCEPTIBILITY AND CAUSALITY 

Since a(u>) is the Fourier transform of a response func- 
tion, it is holomorphic and without zeros in the upper 
complex half-plane (see, e.g., Ref. [2(j for a summary of 
response function properties). Consequently, xm{y,w) as 
given by (the first line of) Eq. i|37fl is there also holomor- 
phic and without zeros, except for possible poles at w- 
values satisfying the equation r)(r)a(oj) / (3eo) = 1. How- 
ever, since a(ui) is the Fourier transform of a response 
function, it is real only on the imaginary frequency axis, 
where it, beginning with a positive value at w = 0, mono- 
tonically decreases with increasing imaginary frequency. 
Hence, if r)(r)a(0)/ (3eo) > 1 is valid, a pole is observed 
and xm(i",w) would fail to be a response function. From 
the requirement that both a(a>) and XM(r,w) be Fourier 
transforms of response functions, it thus follows that 
the condition r](r)a(0) / (3sq) < 1 must be imposed on 
Eq. 127|). 



APPENDIX D: DERIVATION OF EQ. J40l> 

In order to derive Eq. 1)40(1. we return to Eq. I|32() 
together with Eq. (|37fl and recall that, according to 
Eqs. H21f> and l|22l) . the Green tensor Gy(r,r',w) obeys 
the equation 

Gv(r,r',w) = G(r,r',w) 



(B9) 



/ d 3 sG(r,s,w) • x M (s,w)Gv(s,r',w). (Dl) 

JVm 
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According to the principal volume method |3lJ, the de- 
composition 



G(r, s, lo) = VG(r, s,uj) j LS(r — s) 



(D2) 



is applicable to the Green tensor G(r, s,u>) in Eq. (|D1|I . 
The procedure to be followed is to choose first an exclu- 
sion volume whose shape determines the tensor L, and to 
subsequently treat integrals over VG(r,r' ,u>) as (shape- 
dependent) principal value integrals over an exclusion 
volume with the specified shape (see, e.g., Ref. [26j for 
details). Using Eq. (|D2|I in Eq. (|D1|I . and adopting a 
spherical exclusion volume for which L = 1/3, we have 

G v (r, r', iO = [1 + Xu(r, i^)/^ 1 [G(r, r', *£) 
-L f d 3 sVG(r,s,iO-XM(s,iOGv(8,v',iO]- (D3) 

When inserted in Eq. ilML'l) , the second term on the right- 
hand side of Eq. I|D3(I leads to a double integral over r 



and s. Since contributions with s = r are left out in the 
principal value integral, the corresponding part of the 
force is associated with at least two-atom interactions 
in Vm- Dropping all these multi-atom contributions, we 
obtain 



2ttc 2 



/ d 3 r XM (r,;£)VTr 



G(v,v\iQ 
l + XM(r,tO/3 



(D4) 



Applying Eq. I|D4|) to a micro-object whose number den- 
sity of atoms is constant over the small volume Vm, so 
that a position-independent Clausius-Mosotti suscepti- 
bility according to Eq. (|37|l can be assigned to it, and 
replacing [G(r, r', i£)] r '_> r with G( s )(r, r, we just ar- 
rive at Eq. (J301- 
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